Topology of molecular deformations induces triphasic catch bonding in selectin-ligand bonds.

Among the long-standing efforts to elucidate the physical mechanisms of protein-ligand catch bonding, particular attention has been directed at the family of selectin proteins. Selectins exhibit slip, catch-slip, and slip-catch-slip bonding, with minor structural modifications causing major changes in selectins' response to force. How can a single structural mechanism allow interconversion between these various behaviors? We present a unifying theory of selectin-ligand catch bonding, using a structurally motivated free energy landscape to show how the topology of force-induced deformations of the molecular system produces the full range of observed behaviors. We find that the pathway of bond rupture deforms in non-trivial ways, such that unbinding dynamics depend sensitively on force. This implies a severe breakdown of Bell's theory-a paradigmatic theory used widely in catch bond modeling-raising questions about the suitability of Bell's theory in modeling other catch bonds. Our approach can be applied broadly to other protein-ligand systems.

Finite-strain elasticity theory and liquid-liquid phase separation in compressible gels.

The theory of finite-strain elasticity is applied to the phenomenon of cavitation observed in polymer gels following liquid-liquid phase separation of the solvent, which opens a fascinating window on the role of finite-strain elasticity theory in soft materials in general. We show that compressibility effects strongly enhance cavitation in simple materials that obey neo-Hookean elasticity. On the other hand, cavitation phenomena in gels of flexible polymers in a binary solvent that phase separates are surprisingly similar to those of incompressible materials. We find that, as a function of the interfacial energy between the two solvent components, there is a sharp transition between cavitation and classical nucleation and growth. Next, biopolymer gels are characterized by strain hardening and even very low levels of strain hardening turn out to suppress cavitation in polymer gels that obey Flory-Huggins theory in the absence of strain hardening. Our results indicate that cavitation is, in essence, not possible for polymer networks that show strain hardening.

Spanning tree model and the assembly kinetics of RNA viruses.

Single-stranded RNA (ssRNA) viruses self-assemble spontaneously in solutions that contain the viral RNA genome molecules and viral capsid proteins. The self-assembly of empty capsids can be understood on the basis of free energy minimization. However, during the self-assembly of complete viral particles in the cytoplasm of an infected cell, the viral genome molecules must be selected from a large pool of very similar host messenger RNA molecules and it is not known whether this also can be understood by free energy minimization. We address this question using a simple mathematical model, the spanning tree model, that was recently proposed for the assembly of small ssRNA viruses. We present a statistical physics analysis of the properties of this model. RNA selection takes place via a kinetic mechanism that operates during the formation of the nucleation complex and that is related to Hopfield kinetic proofreading.

Packaging contests between viral RNA molecules and kinetic selectivity.

The paper presents a statistical-mechanics model for the kinetic selection of viral RNA molecules by packaging signals during the nucleation stage of the assembly of small RNA viruses. The effects of the RNA secondary structure and folding geometry of the packaging signals on the assembly activation energy barrier are encoded by a pair of characteristics: the wrapping number and the maximum ladder distance. Kinetic selection is found to be optimal when assembly takes place under conditions of supersaturation and also when the concentration ratio of capsid protein and viral RNA concentrations equals the stoichiometric ratio of assembled viral particles. As a function of the height of the activation energy barrier, there is a form of order-disorder transition such that for sufficiently low activation energy barriers, kinetic selectivity is erased by entropic effects associated with the number of assembly pathways.

Physics of viral dynamics.

Viral capsids are often regarded as inert structural units, but in actuality they display fascinating dynamics during different stages of their life cycle. With the advent of single-particle approaches and high-resolution techniques, it is now possible to scrutinize viral dynamics during and after their assembly and during the subsequent development pathway into infectious viruses. In this Review, the focus is on the dynamical properties of viruses, the different physical virology techniques that are being used to study them, and the physical concepts that have been developed to describe viral dynamics.

Finite Temperature Phase Behavior of Viral Capsids as Oriented Particle Shells.

A general phase plot is proposed for discrete particle shells that allows for thermal fluctuations of the shell geometry and of the inter-particle connectivities. The phase plot contains a first-order melting transition, a buckling transition, and a collapse transition and is used to interpret the thermodynamics of microbiological shells.

Invariant theory and orientational phase transitions.

The Landau theory of phase transitions has been productively applied to phase transitions that involve rotational symmetry breaking, such as the transition from an isotropic fluid to a nematic liquid crystal. It even can be applied to the orientational symmetry breaking of simple atomic or molecular clusters that are not true phase transitions. In this paper, we address fundamental problems that arise with the Landau theory when it is applied to rotational symmetry breaking transitions of more complex particle clusters that involve order parameters characterized by larger values of the l index of the dominant spherical harmonic that describes the broken symmetry state. The problems are twofold. First, one may encounter a thermodynamic instability of the expected ground state with respect to states with lower symmetry. A second problem concerns the proliferation of quartic invariants that may or may not be physical. We show that the combination of a geometrical method based on the analysis of the space of invariants, developed by Kim to study symmetry breaking of the Higgs potential, with modern visualization tools provides a resolution to these problems. The approach is applied to the outcome of numerical simulations of particle ordering on a spherical surface and to the ordering of protein shells.

Gaussian curvature and the budding kinetics of enveloped viruses.

The formation of a membrane-enveloped virus starts with the assembly of a curved layer of capsid proteins lining the interior of the plasma membrane (PM) of the host cell. This layer develops into a spherical shell (capsid) enveloped by a lipid-rich membrane. In many cases, the budding process stalls prior to the release of the virus. Recently, Brownian dynamics simulations of a coarse-grained model system reproduced protracted pausing and stalling, which suggests that the origin of pausing/stalling is to be found in the physics of the budding process. Here, we propose that the pausing/stalling observed in the simulations can be understood as a purely kinetic phenomenon associated with the neck geometry. A geometrical potential energy barrier develops during the budding that must be overcome by capsid proteins diffusing along the membrane prior to incorporation into the capsid. The barrier is generated by a conflict between the positive Gauss curvature of the assembling capsid and the negative Gauss curvature of the neck region. A continuum theory description is proposed and is compared with the Brownian simulations of the budding of enveloped viruses.

Specific inter-domain interactions stabilize a compact HIV-1 Gag conformation.

HIV-1 Gag is a large multidomain poly-protein with flexible unstructured linkers connecting its globular subdomains. It is compact when in solution but assumes an extended conformation when assembled within the immature HIV-1 virion. Here, we use molecular dynamics (MD) simulations to quantitatively characterize the intra-domain interactions of HIV-1 Gag. We find that the matrix (MA) domain and the C-terminal subdomain CActd of the CA capsid domain can form a bound state. The bound state, which is held together primarily by interactions between complementary charged and polar residues, stabilizes the compact state of HIV-1 Gag. We calculate the depth of the attractive free energy potential between the MA/ CActd sites and find it to be about three times larger than the dimerization interaction between the CActd domains. Sequence analysis shows high conservation within the newly-found intra-Gag MA/CActd binding site, as well as its spatial proximity to other well known elements of Gag -such as CActd's SP1 helix region, its inositol hexaphosphate (IP6) binding site and major homology region (MHR), as well as the MA trimerization site. Our results point to a high, but yet undetermined, functional significance of the intra-Gag binding site. Recent biophysical experiments that address the binding specificity of Gag are interpreted in the context of the MA/CActd bound state, suggesting an important role in selective packaging of genomic RNA by Gag.

Generalized Flory Theory for Rotational Symmetry Breaking of Complex Macromolecules.

We report on spontaneous rotational symmetry breaking in a minimal model of complex macromolecules with branches and cycles. The transition takes place as the strength of the self-repulsion is increased. At the transition point, the density distribution transforms from isotropic to anisotropic. We analyze this transition using a variational mean-field theory that combines the Gibbs-Bogolyubov-Feynman inequality with the concept of the Laplacian matrix. The density distribution of the broken symmetry state is shown to be determined by the eigenvalues and eigenvectors of this Laplacian matrix. Physically, this reflects the increasing role of the underlying topological structure in determining the density of the macromolecule when repulsive interactions generate internal tension. Eventually, the variational free energy landscape develops a complex structure with multiple competing minima.

Kirigami and the Caspar-Klug construction for viral shells with negative Gauss curvature.

In this work we extend the Caspar-Klug construction to the archaeal viruses, which in recent years have captured the attention of many researchers for their ability to thrive in extreme environments. We assume that the shells of archaeal viruses are composed of hexamers and pentamers-as is true for icosahedral viruses-together with heptamers, necessary to introduce negative Gauss curvature. Following the original work of Caspar and Klug, we first construct models capable of reproducing the shape observed in electron microscopy images of archaeal viruses. Next, using the technique of kirigami, we present a systematic way to formulate archaeal virus templates from regular hexagonal lattices. Finally, we utilize the presented techniques to build finite element models of archaeal virus geometries and investigate their shapes as a function of material properties. In particular, using thin-shell elasticity theory, we describe a buckling transition as a function of a modified Föppl-von Kármán number γ^{★} and we show how changes in γ^{★} may initiate the tail formation in the Acidianus two-tailed archaeal virus.

Amphibian sacculus and the forced Kuramoto model with intrinsic noise and frequency dispersion.

The amphibian sacculus (AS) is an end organ that specializes in the detection of low-frequency auditory and vestibular signals. In this paper, we propose a model for the AS in the form of an array of phase oscillators with long-range coupling, subject to a steady load that suppresses spontaneous oscillations. The array is exposed to significant levels of frequency dispersion and intrinsic noise. We show that such an array can be a sensitive and robust subthreshold detector of low-frequency stimuli, though without significant frequency selectivity. The effects of intrinsic noise and frequency dispersion are contrasted. Intermediate levels of intrinsic noise greatly enhance the sensitivity through stochastic resonance. Frequency dispersion, on the other hand, only degrades detection sensitivity. However, frequency dispersion can play a useful role in terms of the suppression of spontaneous activity. As a model for the AS, the array parameters are such that the system is poised near a saddle-node bifurcation on an invariant circle. However, by a change of array parameters, the same system also can be poised near an emergent Andronov-Hopf bifurcation and thereby function as a frequency-selective detector.

Structural studies of Acidianus tailed spindle virus reveal a structural paradigm used in the assembly of spindle-shaped viruses.

The spindle-shaped virion morphology is common among archaeal viruses, where it is a defining characteristic of many viral families. However, structural heterogeneity intrinsic to spindle-shaped viruses has seriously hindered efforts to elucidate the molecular architecture of these lemon-shaped capsids. We have utilized a combination of cryo-electron microscopy and X-ray crystallography to study Acidianus tailed spindle virus (ATSV). These studies reveal the architectural principles that underlie assembly of a spindle-shaped virus. Cryo-electron tomography shows a smooth transition from the spindle-shaped capsid into the tubular-shaped tail and allows low-resolution structural modeling of individual virions. Remarkably, higher-dose 2D micrographs reveal a helical surface lattice in the spindle-shaped capsid. Consistent with this, crystallographic studies of the major capsid protein reveal a decorated four-helix bundle that packs within the crystal to form a four-start helical assembly with structural similarity to the tube-shaped tail structure of ATSV and other tailed, spindle-shaped viruses. Combined, this suggests that the spindle-shaped morphology of the ATSV capsid is formed by a multistart helical assembly with a smoothly varying radius and allows construction of a pseudoatomic model for the lemon-shaped capsid that extends into a tubular tail. The potential advantages that this novel architecture conveys to the life cycle of spindle-shaped viruses, including a role in DNA ejection, are discussed.

Confining annealed branched polymers inside spherical capsids.

The Lifshitz equation for the confinement of a linear polymer in a spherical cavity of radius R has the form of the Schrödinger equation for a quantum particle trapped in a potential well with flat bottom and infinite walls at radius R. We show that the Lifshitz equation of a confined annealed branched polymer has the form of the Schrödinger equation for a quantum harmonic oscillator. The resulting confinement energy has a 1/R4 dependence on the confinement radius R, in contrast to the case of confined linear polymers, which have a 1/R2 dependence. We discuss the application of this result to the problem of the confinement of single-stranded RNA molecules inside spherical capsids.

Ground state instabilities of protein shells are eliminated by buckling.

We propose a hybrid discrete-continuum model to study the ground state of protein shells. The model allows for shape transformation of the shell and buckling transitions as well as the competition between states with different symmetries that characterize discrete particle models with radial pair potentials. Our main results are as follows. For large Föppl-von Kármán (FvK) numbers the shells have stable isometric ground states. As the FvK number is reduced, shells undergo a buckling transition resembling that of thin-shell elasticity theory. When the width of the pair potential is reduced below a critical value, then buckling coincides with the onset of structural instability triggered by over-stretched pair potentials. Chiral shells are found to be more prone to structural instability than achiral shells. It is argued that the well-width appropriate for protein shells lies below the structural instability threshold. This means that the self-assembly of protein shells with a well-defined, stable structure is possible only if the bending energy of the shell is sufficiently low so that the FvK number of the assembled shell is above the buckling threshold.

Orientational phase transitions and the assembly of viral capsids.

We present a Landau theory for large-l orientational phase transitions and apply it to the assembly of icosahedral viral capsids. The theory predicts two distinct types of ordering transitions. Transitions dominated by the l=6,10,12, and 18 icosahedral spherical harmonics resemble robust first-order phase transitions that are not significantly affected by chirality. The remaining transitions depend essentially on including mixed l states denoted as l=15+16 corresponding to a mixture of l=15 and l=16 spherical harmonics. The l=15+16 transition is either continuous or weakly first-order and it is strongly influenced by chirality, which suppresses spontaneous chiral symmetry breaking. The icosahedral state is in close competition with states that have tetrahedral, D_{5}, and octahedral symmetries. We present a group-theoretic method to analyze the competition between the different symmetries. The theory is applied to a variety of viral shells.

Sequence Dependence of Viral RNA Encapsidation.

We develop a Flory mean-field theory for viral RNA (vRNA) molecules that extends the current RNA folding algorithms to include interactions between different sections of the secondary structure. The theory is applied to sequence-selective vRNA encapsidation. The dependence on sequence enters through a single parameter: the largest eigenvalue of the Kramers matrix of the branched polymer obtained by coarse graining the secondary structure. Differences between the work of encapsidation of vRNA molecules and of randomized isomers are found to be in the range of 20 kBT, more than sufficient to provide a strong bias in favor of vRNA encapsidation. The method is applied to a packaging competition experiment where large vRNA molecules compete for encapsidation with two smaller RNA species that together have the same nucleotide sequence as the large molecule. We encounter a substantial, generic free energy bias, that also is of the order of 20 kBT, in favor of encapsidating the single large RNA molecule. The bias is mainly the consequence of the fact that dividing up a large vRNA molecule involves the release of stored elastic energy. This provides an important, nonspecific mechanism for preferential encapsidation of single larger vRNA molecules over multiple smaller mRNA molecules with the same total number of nucleotides. The result is also consistent with recent RNA packaging competition experiments by Comas-Garcia et al.1 Finally, the Flory method leads to the result that when two RNA molecules are copackaged, they are expected to remain segregated inside the capsid.

Impurities in Bose-Einstein Condensates: From Polaron to Soliton.

We propose that impurities in a Bose-Einstein condensate which is coupled to a transversely laser-pumped multimode cavity form an experimentally accessible and analytically tractable model system for the study of impurities solvated in correlated liquids and the breakdown of linear-response theory [corrected]. As the strength of the coupling constant between the impurity and the Bose-Einstein condensate is increased, which is possible through Feshbach resonance methods, the impurity passes from a large to a small polaron state, and then to an impurity-soliton state. This last transition marks the breakdown of linear-response theory.

Phase-locked spiking of inner ear hair cells and the driven noisy Adler equation.

The inner ear constitutes a remarkably sensitive mechanical detector. This detection occurs in a noisy and highly viscous environment, as the sensory cells-the hair cells-are immersed in a fluid-filled compartment and operate at room or higher temperatures. We model the active motility of hair cell bundles of the vestibular system with the Adler equation, which describes the phase degree of freedom of bundle motion. We explore both analytically and numerically the response of the system to external signals, in the presence of white noise. The theoretical model predicts that hair bundles poised in the quiescent regime can exhibit sporadic spikes-sudden excursions in the position of the bundle. In this spiking regime, the system exhibits stochastic resonance, with the spiking rate peaking at an optimal level of noise. Upon the application of a very weak signal, the spikes occur at a preferential phase of the stimulus cycle. We compare the theoretical predictions of our model to experimental measurements obtained in vitro from individual hair cells. Finally, we show that an array of uncoupled hair cells could provide a sensitive detector that encodes the frequency of the applied signal.